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ABSTRACT 



A study of an Automatic Power Compensator with inputs of airspeed 
and angle of attack was conducted by means of a simplified analytical 
analysis and analogue simulation. Digital Computer programming was 
used for data reduction. The Chance Vought F-8 aircraft was used as 
the vehicle for the study. Loop gains were initially set by means of a 

I# % 

simplified analysis based on the independent phugoid and short period 
motions. The validity of such an analysis was demonstrated by ana- 
logue simulation of the complete system. The system is feasible^.bdth 
for independent use or as a component of a complete automatic flight 
system. 
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TABLE OF SYMBOLS 



(For those symbols not defined in the text) 
ac Aerodynamic center 

B (Iyy) foment of Inertia about the Y axis 

c chord (MAC of wing) 

Cl coefficient of Lift 

Cj) coefficient of Drag 

C m coefficient of Moment 

eg center of gravity 

g acceleration of gravity 

K* angle of attack proportional gain 

K u airspeed proportional gain 

k u K u airspeed integral gain 

L Lift force 

1^, j? t tail length - distance from eg to ac of horizontal tail 
m mass 

n vertical acceleration perturbation (g units) 

u horizontal velocity perturbation 

U D steady state velocity 

V reference velocity (V = U Q ) 

w vertical velocity perturbation 

W aircraft weight 

S w reference area (wing area) 
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s 

T 

(t) 

oc 

V 

e 

7 

P 

T 



Laplace operator 
thrust 

indicates variable in time domain 
angle of attack, (see Fig 2) 

U 0 

elevator deflection (see Fig 2) 

angle above horizontal (see Fig 2) 

flight path angle (see Fig 2) 

atmospheric density 

time constant for delays (see Table IV) 

indicates determinant expansion 



indicates matrix formulation 



Superscripts: 

-d i I 

dt 



dii 1 

d t 2 



Subscripts: 

o initial or steady state value 



e 

c 

f 



error quantity or signal 
command signal 
delayed signal 



Aerodynamic derivatives 

The standard NACA notation is used for the non-dimensional derivatives 
in that v ^ 

£> c _l. c _ oc, 

< '-'I 



CLu = 



6 U 



'm 



m 



oc 






; etc. (See Table II) 



12 



1 . 



Introduction. 



This study is a continuation of a series of investigations designed 
to evaluate various types of "Auto -throttle" devices. The study was in- 
stituted in the Aeronautical Engineering Department of the United States 
Naval Postgraduate School in 1963. (8, 9) 

The need for such a device on high performance carrier aircraft and 
the basic concepts relative to the system have been previously discussed. 
Briefly, the power compensator is a device meant to automatically stabi- 
lize the airspeed at a pre-selected value without attention from the pilot. 
The pilot then maintains the desired attitude or glide slope. This concept 
is not limited to, but is especially desirable during carrier type approaches 
where modem, high performance, swept-wing aircraft must operate in a 
range of inherent speed instability. With the heavier aircraft now being 
programmed for carrier use lower landing speeds are most dppirable. By 
allowing the pilot to focus his attention on only one correction during ap- 
proach, . a power compensator permits the safe utilization of slower approach 
speeds. Basically then, the automatic power compensation is a means of 
artificially tailoring the aircrafts inherent flight characteristics in the low 
speed region to meet the opposing requirements of safety and lower landing 
speeds. In this sense such a device could prove applicable, even neces- 
sary, to the landing of the larger super-sonic transport aircraft now under 
development. 

The initial study demonstrated the need for such a device by showing 
real time, pilot controlled, analogue simulated mirror approaches with and 
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without an auto-throttle. [ 8 ] The second paper attempted to improve the 
systems by using different inputs. [9] The systems analyzed thus far use 
inputs of airspeed/pitch angle and angle of attack/normal acceleration. 
These are systems of the type investigated by the United States Navy for 
operational use. The angle of attack/normal acceleration system is pre- 
sently in use on operational aircraft. 

Utilizing the results of the previous studies. Professor E. J. Andrews 
of the U. S. Naval Postgraduate School proposed a much simplified analytic 
means of setting loop gains. [4] 

This paper will investigate a power compensator that uses inputs of 
airspeed and angle of attack. Within the knowledge of the author no such 
system is in use at present. This type of a system might be thought of as 
a "hybrid system" since it utilizes one of the inputs from each of the two 
systems investigated thus far. The motivation for such a system lies in the 
fact that airspeed and angle of attack are two of the most "basic" governing 
parameters of flight. Since these parameters are related by 
( 4 + m)W = ^pV 2 SC LoC <=C 

4U 

the values of airspeed and angle of attack will be related to the aircraft 
weight and the desired type of flight path. For a given aircraft weight, the 
longitudinal .flight. path is determined completely by the elevator position and 
the throttle setting (airspeed). 

Therefore, with such an auto-throttle stablizing airspeed, the pilot 
(or autopilot) need only control the pitch attitude to obtain any desired flight 
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path. This type of an auto-throttle would be especially useful for flying 
a pre- determined flight path designed to give maximum range, maximum 
rate of climb or a desired glide path. over a range of airspeeds. 

In the investigation of this type of a power compensator, the simpli- 
fied analytic approach will be demonstrated in approximating the desired 
gain settings. The Chance-Vought F-8 aircraft was. used as the vehicle in 
the previous studies and will be used in this study. 

The equations of motion are solved using an analogue computer. A 
digital computer is used to determine aircraft parameters and roots of char- 
acteristic equations. 

All aerodynamic parameters are used in the English dimensional form. 
All equations and programming are given in real time. 

The aerodynamic data used was obtained from the Chance Vought 
Division of Ling-Temco Vought and an earlier report on an evaluation of a 
power compensator. [6, 7] 

A word of appreciation is due Professor E. J. Andrews for his assis- 
tance during this investigation. 

2. Discussion. 

2. 1 General 

The system under study is actually equivalent to the angle of attack/ 
normal acceleration system evaluated previously. The auto-throttle equa- 
tion for the oc/n system is: 

Eq(2 — l). T c = K^oc — I^n 
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Since 



L - -lp(Vit:jl» 2 SC L = ±-p. SC L (V 2 i* 2uV3 u2) 

Then = -|pv 2 SC LlC ; = ^pSC L (2V+ 2u) ^ SC L Vp 

nW = AL 

6 n _ 1 & L . & n _ _1_ _^_L 

Oo c W £)<*: * 6 u W ^ u 

Therefore: 

n = _ ^P v2 ~w C l «- cC + P “v7 C L Vu = £ 3 °° " K 4 u 

The Auto-throttle equation becomes: 

T c = K : oC - K 2 (K 3 oC + K 4 u) = (K l - K 2 K 3 )oC- K^u 

Eq(2-2) T c = K^oC - K u u 

A block diagram of the complete system is shown in Figure 1. 

2. 2 Airframe Equations 

The 3 -degree of freedom, small perturbation, aircraft longitudinal 
equations were used to solve for the airframe motion. The equations shown 
below are as developed in any standard text on aerodynamic stability. [3] 
The definitions of the derivatives are listed in Table I. The dimensional 
derivatives were computed from the non-dimensional derivatives utilizing 
a digital computer program given in the Appendix. The equations used for 
this conversion are given in Table II. The values of the dimensional deri- 
vatives obtained are listed in Table III along with all other pertinent flight 
parameters. 

Eqs(2-3) u = X u u + X w w + X q 0 - gcos0 o 0 + X T T + X^T) 

(l - Z^)w = Z y u + Z w w + (Z q + U 0 )0 - gsin0 o 0 + Z^T + Z tf~) 
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0 = M u u + M w w + M^w + Mq0 + M>pT + M rfTj 
The angular and velocity relationships are defined in Figure 2. The 
derivatives concerning thrust are dervied in Figure 3. These equations were 
reduced by inserting the values ; 




Z- 



w 



0 ; 



oC = 



w 

Uo 



and rearranged to give 

Eqs(2-4) u = X u u + U Q X W «< - gcos0<,0 + X-pT + 

*--§tu + z w < + e--a^-e + -^T + ^j 

0 = M u u + U q M w «< + U q M^cC + M q 0 + M t T + UyfTj 



2.3 Power Compensator 
The equation of the power compensator is: 



Eqs(2-5) 



T„ = 



1 + % s 






K 



JJL 



1 + 



■(1 + 






■)u. 



where % and % are estimates of the angle of attack and airspeed measure- 
ment delays respectively. It can be seen from this equation that the power 
compensator will apply a forward ( + ) thrust for a nose up (+) angle of attack 
perturbation and an increase in airspeed (+ u e ) will cause a reduction of 
thrust. The compensator utilizes angle of attack and airspeed proportional 
feedbacks plus an airspeed integration to insure zero steady state errors in 
the airspeed loop. 
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2. 4 Transfer Functions 



The transfer functions used are defined in Table IV. The delays were 
obtained from various sources. [6,7] The airframe transfer functions are 
derived in the Appendix from the airframe Laplace equations shown below 
in matrix form. The airframe transfer functions were evaluated by means 
of a digital computer program given in the Appendix* 



Eqs(2-6) Is - X u 



- U o X w 



- Z il 



- M 



yt 



+ (s - z w ) 

- U 0 (m w + sm^) 



+ gcos0 Q 

_ s + gsin 3 Q 

U 



\ 



u \ 



S2 - MS I ' 0 j 





Xij 


X * \ 




= 


Zx 


_2a_ 






Uo 


Uo 






m t 


M r\ / 


\ 



V 



The open and closed loop transfer functions were derived from 
Figure 1 using standard methods. 

V ‘A* 

Open Loop: .. ' '' <' „• 



Eqs(2-7) 



u _ 



Ur 



D- T ° U 

- D U rr- - 

u ci 1 



-±£_ 






! + 1L 1 a Is ± 

T c ^ ^ * T 



v 



Eqs(2-8) 



u _ 



|c |f T_ ^ 
u e T g Tf T 



u 



Ucb Tfc Tf T ^ T fc Tf T 



T c Tf x <<e 



Only the airspeed transfer functions are given above. Transfer func- 
tions for ^/a Ce . , u /tj » <*•/ '-q could also be developed easily from the above 
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expression since the denominator, ie. the characteristic equation, would 
be the same in all cases. 

2. 5 Digital Computer Programs 

The Digital Computer programs used in the previous studies were 
rewritten in FORTRAN 63 in a subroutine format. These programs are used 
to convert either wind tunnel data or non-dimensional derivatives to the 
dimensional derivatives necessary for the Analogue Program. These deriv- 
atives are then used as inputs to other subprograms which analize the basic 
airframe motion and cpmpute the airframe transfer functions for the inputs 
of elevator and thrust. Thus it is possible to proceed directly from wind 
tunnel data to the response of the airframe for several flight conditions, 
using these results to analize any gain changes necessary in the control 
loops. These programs are contained in the Appendix. 

2.6 Complete Analysis 

The above transfer functions, both open and closed loop, could now 
be used to analyze the system. This involves a 5th order characteristic 
equation with at least three variable gains. The numerators are of third 
order or less. The characteristic equation is evaluated in the Appendix 
for later use in a Root Locus Diagram. At this point the methods of Root 
Locus plots combined with frequency response methods could be used to 
analyze the response for various combinations of gains to elevator inputs, 
and step gusts of u and oC. This procedure would require an attack by 
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digital computor methods for at least computing the roots of the various 
polynomials generated for each value of the gains used and then perhaps 
estimating the response obtained by means of simplified Bode plots or 
Root Locus plots. But this rigorous, automatic control engineering appro- 
ach at this point is unnecessary and uneconomical if not also mathematic- 
ally unsound. Within the accuracy of the assumption inherent in the basic 
airframe equations (Eqs. 2-3) a much simplified approach is justified. In 
fact, it is known that the responses calculated by this "exact" method are 
not only in error, but for the case of shear gusts, can give un conservative 
results for structural loading parameters. [3] 



2.7 Simplified Analysis 

This analysis is based on the standard assumption of reducing the 
airframe equations of motion by assuming two independent motions; one at 
constant angle of attack, and the other at constant airspeed. Usintj.this ajssunp- 
tion, Eqs. (2—6) were separated and reduced as shown below: (negligible 
values from Table III are omitted) 

For Phugoid motion with cC= — S / n ^ Q =0 ; cos0_ = 1.0 

J ..Uo ° 



Eqs(2-9) 



s - x u 


g 




u 




x T 


0 1 




' T 


- Zu 


- s 




e 




vO 


z* 




V 


u 0 












'Uo , 
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Then 



Eqs(2-10) 



Xj 
0 



g 

- s 



SXn 



-z, 

u. 



'u 



u 



- - s- 



S 2 - X..S - g -5lL- 
u y U, 



Eqs(2-ll) 



- g 



U 

v 



u. 



S 2 - x u s - 



2u 



For the Short Period Motion with u = 0 
Eqs(2-12) /S - Z w - S 



cc 



i- U Q (M w + SM^) S 2 - M q S / \ 0 



Z*. 

U, 



M 



T ^*1 / \ ^7 



Eqs(2-13) 



0 

M^ 



- S 



S2 - 



MqS 



S - Z. 



w 



- S 



-U Q (M w + SM W ) S 2 - 



M q S 



Eqs(2-14) 



Mn 



s 2 - (Zw + Mq + U 0 M*)S+Z w M q - U Q M W 



Mi 



MqZ^ 

U„ 



V "si - (Z w +' M q + U 0 M^ ) S + Z w M q - U 0 M W 



In this way the aircraft equations of motion reduce to second order 
Laplace equations which are extremely simple to handle. The time response 
to any type of input of either elevator or thrust can now be approximated by 
standard Laplace techniques. The response of airspeed to either thrust or 
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elevator will be controlled by the lightly damped (at this airspeed) phugoid 
characteristic equation while the angle of attack, response to either forc- 
ing function will have the characteristic of the short period mode. 

At this point a further simplification can be made by examining the 
short period motion. By substituting the aerodynamic values in the above 
equations it can be seen that the short period mode has a damping ratio of 
. 35 with a time to one tenth of the initial disturbance .vaiue obcSu.&seteonds 
tre&qieriod o£<£bfleconds. 

Thus it can be seen that the response to this disturbance will, ap- 
•preaohchnuexpQneintially damped time response rather than an oscillatory 
type. This leads directly to the well known assumption that the aircraft 
moment of inertia in the pitch plane can be neglected. In the Laplace 
equations above it can be reflected by eliminating the S2 term from the 



matrix. 



1 



The equations then reduce to: 
Eqs(2-15) 



M „ -512k. 



TT 



o C ^ 

T ~ S ( - U 0 M^, - Mq) + Z^M, - U C M W 

oc _ 



I S ( - U 0 M* - M q ) + Z w M q - U q M, 

Which are of the form 

oC C 9 



’1 



Eqs(2-16) = __ 

7] S + A 



S + A 



1, 



This is more evident from the NACA nondimensional form of the equa- 
tions. The coefficient of the term is then much smaller than the other 
terms. 
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where A is 



A = 

-U q M^ - M q 



For a step input of elevator or thrust, Eqs(2-16) transform to the time do- 
main as 

Eqs(2rl7) °C(t) = £i(. ^~ At \ ^ (t) _ £l2 

r){t) A \ ‘ ° j T(t) A 

Since A is usually large (3.5 in this case) the exponentionai term can be 

neglected. This indicates that the angle of attack response to either thrust 

or elevator will be approximately constant. In aerodynamic terms, 

Eqs(2-18) ^(t) M » - Z, M g /U 0 oC( t ) _ M T 

Tjb.) Z w M q - U 0 M W T(t) Z w M q - U 0 M W 




Substituting values gives, 
Eqs(2-19) £ (t) _ 

V (t) " " 



= 

T(t) 



-1.2 x 10 ^rad/lb 



The value obtained for oC(t)/TJ (t) is within 10% of the value deter- 
mined previously by analogue results. [9] The low value for <^c(t)/T(t) shows 
that the angle of attack response to thrust is negligible. The Power Com- 
pensator then can be used only to control the airspeed perturbations ^and 
elevator must be used to control the angle of attack perturbations. This is 
not to say that the airspeed variations will not affect the angle of attack. 
Any airspeed perturbations will cause the aircraft to seek a new equilibrium 
angle of attack. The angle of attack loop is necessary to offset the large 
rise in drag for positive excursions of angle of attack. 
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The assumption of two independent motions is now carried into the 
control system of Figure jL by considering each loop separately. The equa- 
tions are further simplified by assuming all delay terms are equal to unity. 
The two simplified loops with their closed loop transfer functions are shown 
in Figure 4. These functions with aerodynamic terms substituted are shown 
below. 

Eqs(2-20) _u = * k u*. 

u c S2 + (K u X t - 3^)8 - Z u + XTK u k u 

Eqs(2-21) = ifl 

oC c S2 + (Z w + M q )S + Z Vf M q - U q M w - 

Roots Locus plots for these functions can be easily estimated since 
they are only second order. Root Locus plots for these functions are shown 
in Figures 5 and 6 for the aerodynamic values listed in Table III. 

The gains can then be chosen from these Root Locus Plots. As a first 
approximation the gains can be chosen such that the phugoid roots are on 
the imaginary axis in close proximity to the short period roots. This will 
give a response of moderate overshoot with the damping of the short period 
roots. Note that the short period roots are changed only slightly from their 
basic aircraft position and that the gain of the angle of attack loop affects 
these roots only to a small degree. 

The angle of attack gain can be approximated by considering its de- 
sired function of equalizing the drag rise due to the angle of attack. This 
gives 
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where To = K^oC-e 



b T c = b D 

boC b<*~ 

Therefore K,c = — pV^SCj-^. This gives an approximation of 24, 000 lb/ 

2 

rad. The values of the other gains chosen were : 

K u = 300 - 500 lb/ft/sec 
k u = . 05/sec 

The effect of the gains on the overall motion can be analized from 
the simplified Root Locus. K u will mainly affect the time constant of the 
response. Thus increasing K^will increase the speed of the response with 
a resulting increase in overshoot. The effect of K«c is small but will increase 
the overshoot as it is increased. The effect of k u is difficult to analize 
without a detailed study of the Root Locus. In general, increasing k u will 
increase the overshoot and can cause a second oscillatory mode to appear 
at about the same frequency as the short period mode. This would give an 
undesirables response. 

The gains chosen as first approximations are : 

K u = 400 lb/ft/sec 
k u = . 05/sec 

K<. = 20, 000 lb/rad 
K u k u = 2.0 lb/ft 

2.8 Co-ordinated Elevator Input 

In order to stabalize the aircraft motion at a steady value of airspeed, 
elevator must be used to control the flight path. This input will simulate 
the action of either an autopilot feedback or the action of the pilot. The 
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gain for this elevator/airspeed loop is obtained from the following analysis; 
V2 = (U. + u)2 = = U 2 + 2uU 

ipct 



then 



or 



2U du 



- W/S 




dC L 



2uU 



- W/S 

Jp c i 



C L^ oC 



therefore 



u 



-W/SC Tjt ^ 
PUC L 2 



For the given flight condition 

u(kts) = -140. 0 °£(rad) = -6.32oC(deg) 
From Eqs 19, = — 1.74 ^7 

Then u(kts) = 11. 0 7 ? (deg) 



3. Analogue Simulation. 

In order to check the validity of the foregoing analysis, the full 
equations (Eqs. 2-4) were solved on an analogue computer and the gains 
of the Power Compensator "optimized" for horizontal and vertical gust in- 
puts. A Donner 3100 computer was used in the analysis. The analogue 
wiring diagrams are shown in Figures 7 and 8. All analogue program infor- 
mation is given in the Appendix. The values given for the flight condition 
in Table III were used in the program. Figures 9, 10, and 11 show the 
response of the basic airframe to horizontal and vertical gusts and to a 
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step elevator input. The motion shows the phugoid response with a period 
of 28 seconds which compares with the calculated value of 30 seconds. The 
short period motion is difficult to see because of the heavy damping but the 
period is less thpn 10 seconds. Figure 11 clearly shows the relation be- 
tween elevator input and angle of attack. The angle of attack rapidly takes 
on a small phugoid oscillation about a constant value of - 1« 8 deg for an 
elevator input of 1 deg. This compares well with Eq(2-19). 

3.1 Basic Airframe 

Figures 12, 13, and 14 show the response of the basic airframe to 
sinusoidal horizontal gusts with periods close to those of the short period 
and phugoid motion. Figures 12 and 13 show that the response to the short 
period gust is of small magnitude and composed of the short period motion 
superimposed on the longer period phugoid. As the period of the phugoid 
mode is reached. Figure 14 shows that the response becomes divergent. The 
same is true for vertical gust inputs shown in Figures 15 and 16. 

3. 2 Basic Airframe with Power Compensator 
In order to check the validity of the gain analysis given in 2.7, the 
response of the airspeed motion for various values of loop gains is given in 
Figures 17 through 20. Figure 17 shows that K u affects the damping of the 
motion and is the main controlling gain. Figure 18 -shows the decrease of 
damping with increasing integral gain, K u k u . Figures 19 and 20 show the 
small effect of the angle of attack gain, K ^ , on both horizontal and vertical 
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gusts. Since the power compensator will have the most effect on the air- 
speed loop, the airspeed response was used to choose the "best” gains. 

Figures 21 through 25 show the response of all the important para- 
meters to various gains for the most unfavorable gust condition; - 5 kt 
horizontal and + 5 kt vertical gust. Important aspects of these Figures are 
summarized below. 







TABLE V 
ANALOGUE DATA 




Fiq 


Overshoot 


Time to 
within 1 kt 


Maximum T 
- reauired 


'Y after 
10 sec 


21 


65% 


6. 6 sec 


2000 lbs 


-1.5 deg 


22 


80 


8.5 


1800 


- 1.3 


23 


70 


4.8 


2800 


- 1. 5 


24 


85 


9. 5 


3500 


-1.5 


25 


80 


9. 0 


1400 


- 1.2 



The gains from Figure 21 were considered as the most favorable. The over- 
shoot is somewhat misleading since it will be reduced considerably by the 
pilots application of elevator to correct the flight errors. However, this 
overshoot could be reduced by series compensation. 

The gains of Figure 21 were used for all the remaining runs containing 
the power compensator and are listed below. 

K u = 400 Ib/ft/sec 
K u k u = 2.0 lb/ft 
K«c = 10,000 Ib/rad 

Figures 26 through 30 show the response of the airframe and power compen- 
sator to sinusoidal gust inputs. They show that the compensator eliminates 
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the phugoid divergence but Figure 29 indicates that the system tends to 
magnify the response for vertical gusts with a period of 9. 6 seconds. How- 
ever this period should be long enough to be controlled by the pilot and the 
airspeed perturbation is not serious. Figures 31 and 32 show the response 
to an airspeed command with coordinated elevator input as discussed in 2. 8. 
These Figures show that the overshoot is reduced to 26% and the glide slope 
variation to within . 5 .degrees. 

4. Root Locus Analysis . 

In order to compare the simplified Root Locus plots with those for the 
complete system, the characteristic equation for the complete system was 
evaluated from Eq(2-8). This is shown in the Appendix. Root Locus plots 
of this polynomial were then developed for several values of loop gains. 
These plots are shown in Figures 33 and 34. 

5. Conclusions. 

1. The airspeed/angle of attack power compensator would be a feas- 
ible system to employ especially where the use of a full auto-pilot is anti- 
cipated. 

2. The simplified method of setting loop gains by considering the 
control loops as separate functions of the phugoid and short period motions 
will give valid results. 

3. The response of the system can be approximated by considering 
the separate root locus plots of each loop. 

4. The system may have to be improved because of the excessive 
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overshoot response of airspeed. This could be done by series compensation. 

* 

5. The effect of the gains on the Root Locus of the complete system 
compares qualitatively to the effect of the gains on the two loops consider- 
ed separately. 

6. The power compensator will have negligible effect on the short 
period motion and can be thought of as an airspeed (phugoid) damping device. 
This will be true whenever the thrust axis is close to the body axis giving 
small values for M-j, and Zj. 

6. Recommendations for Future Study. 

1. The effects of the non-linear aspects of on the pitch motion 

feedback gain requirements for all three power compensators studied thus 
far.should be the next field of study. 

2. An investigation to determine which of the three systems studied 
thus far would give the best results when combined with an elevator auto- 
pilot would indicate an optimum combination for a complete system. 

3. A complete analysis of a flight control system should include the 
effects of turbulent air gusts and aeroelasticity. To do this type of an analy- 
sis by analogue simulation would require an excessive number of compon- 
ents. Since the problem of response to random frequencies is best analyzed 
by frequency response methods, a digital computer program would give the 
most convenient approach. This type of a problem could also be analyzed 
on a "hybrid computer system" by using the analogue to solve the equations 
and the digital to introduce the non-linearities. 
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TABLE I * 

f 

DEFINITIONS OF DIMENSIONAL AERODYNAMIC DERIVATIVES 



X 




_1_ 




1 


A u 




m 


du 


sec 


Xw 




1 


M. 

c)W 


1 




m 


sec 






1 


dX 


ft 


x q 




m 


’ dq 


sec-rad 


Xp 




1 


DX 


ft 




m 


' ST" 


sec 2 -lb 






1 


dX 


ft 






m 


d>? 


sec 2 -rad 
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dZ 
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Zu 
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du 


sec 


Zw 




, 1 


az 


1 
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t)w 


sec 


Z q 




_1_ 


dZ 


1 




m 


dq 


sec-rad 


Zip 




1 


az 


ft 
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*St* 


sec 2 - lb 


Zn 






az 


ft 




m 




sec 2 -rad 


M u 
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dM 


1 




B 


du 


sec-ft 


Mw 
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dM 


1 
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dw 


sec-ft 
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dM 
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TT 


dw 


ft 


M q 
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dM 


1 
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dq 


sec-rad 


Mp 
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dM 

ST 
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B 


sec 2 -lb 
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dM 
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sec 2 - rad 
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TABLE XI 

CONVERSION OF DERIVATIVES 



Xu 


. 


^“ Cd “ Cdu ^ 


Xv* 


= 






= 




Xj 


= 


COS^/ 

m 


Zu 


= 




Zv, 


— 


(-c L * - Cp) 


Z >1 


= 


^(-cl,) 


Zy 


— 


sinV 

m 




= 


^(-C T z t ) 


% 


= 


*W- to— > 


«* 


= 


«W ^ 


M q 


= 






= 


>^SV2c , x 

2B v ^ 1 



TABLE III 



AIRCRAFT PARAMETERS 

s w 375 ft 2 

ziji — .4- 37 ft 

l t 14 ..O 8 ft 

c 11.78 ft 

B 96,000 slug-ft 2 

FLIGHT CONDITIONS 

V 1 . 15 V S 234 - ft/sec 

eg 2k.fo MAC 

V/ 22,000 lbs 

6 0 8.1 deg 

AERODYNAMIC DERIVATIVES 
-.060 l/sec 
Xv, -. 014-19 l/sec 

X q 0 

X T . 0014-5 ft/lb-sec 2 

X-^ -1. 64 ft/sec 2 -rad 

-.2655 l/sec 
-#[|.265 l/sec 
0 

Z T - 2.17 x 10-5 ft/ sec 2 -lb 

Z>j - 19. 1 ft/sec 2 -rad 

Mu .000185 l/sec-ft 

My, — .004858 l/sec-ft 

-1.772 x 10"^ 1/ft 
Mq -.3384- 1/sec-rad 



3 * 



TABLE III (CONT) 



. i 

'X y 

M t x 10“ 6 l/lb-sec 

M>j -2*25> l/rad-sec 2 



U 
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TABLE IV 



TRANSFER FUNCTIONS 



Airframe 



Characteristic equation (CH) 



sk + .8659S 3 + 1.3.3S 2 + .0608s + .0423 

u _ .00l45(S 3 + .80£9S 2 + 1.378S + .020) 
T “ CH 

u_ .198(S 2 + ljli3.9S + 161.68) 



oc ~6.i6xio~^(s 2 + .0772s + .0267) 

T CH 

cC -.05q6(S 3 + lg.93S 2 + 1.70S + 78.23) 



Engine 

_ 2 _ - 1 

T f 1 +T« S = r l.l£ 

Angle of Attack Measurement Lag 



= 1 S 



T^= 0.5 sec 



Airspeed Measurement Lag 



1 + T 6 s 



T s ~ 0.1 sec 



Power Compensator 







Power Compensator Servo Lag 
Tf 1 



7f = 0.2 sec 



T c ” 1 + T f S 
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BLOCK DIAGRAM OP COMPLETE SYSTEM 





FIGURE 2 



VELOCITY AND ANGULAR RELATIONSHIPS 
Arrows show positive direction 
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FIGURE 3 

DERIVATION OF THRUST DERIVATIVES 
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FIGURE 4 



SIMPLIFIED LOOPS 
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FIGURE 5 



ROOT LOCUS 

Simplified Airspeed Loop - Varies 
Basic Airframe Phugoid Roots A 




IMAGINARY 



FIGURE 6 



ROOT LOCUS 

Simplified Angle of Attack Loop - K* Varies 
Basic Airframe Short Period Roots A 
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ANALOGUE SIMULATION CIRCUITS - Airframe 



ANALOGUE SIMULATION CIRCUITS 
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Basic Airframe 




FIGURE 10 

Basic Airframe -5 kts«=2.1°) Vertical Gust Input 
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FIGURE 12 



Basic Airframe - Sinusoidal Horizontal Gust Input 

P = 5.6 sec 
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Basic Airframe - Sinusoidal Horizontal Gust Input 

? =20 sec 
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Basio Airframe - Sinusoidal Vertical Gust Input 
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CLEVITE corpo 



NATION 



CLEVELAND. OHIO 



PRINTED IN U.S A. 




Airframe and Power Compensator -5 lets Horizontal 
and +5 kts(<=2*l° ) Vertical Gust Inputs 
K^lj.00 1^=2 K*= 1^0,000 
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and +5 kts(oC=2,l°) Vertical Gust Inputs 

1^=600 =.$ Kot — 1|0,000 
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Airframe and Power Compensator -5 kts Horizontal 
and +5 kts(oc=2.1°) Vertical Gust Inputs 
=200 ^ 1 ^ = 2 K* = 10,000 



6l 




Airframe and Power Compensator 
Sinusoidal Vertical Gust Input P = 9*6 sec 
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FIGURE 33 



ROOT LOCUS 

Complete Characteristic Equation 
= 600 K*. = 10,000 k u Varies 




X -SCALE " 2.00E-O1 UMITS^INCH. 
Y -SCALE - 2.00E-01 UNITS'INCH. 



FIGURE 34- 



ROOT LOCUS 

Complete Characteristic Equation 
Ku = lj.00 K* - 10,000 k u Varies 




X -SCALE = 2.09E-01 UNITS/INCH. 
Y -SCALE = 2.00E-01 UNITS'INCH. 
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7. 



APPENDIX 



7. 1 Digital Computer Programs 

The digital computer programs are designed to compute the character- 
istic equation and the transfer functions for elevator and thrust inputs. By 
utilizing the proper combination of subprograms either wind tunnel data and 
airframe dimensions or reported nondimensional derivatives may be used in 
these programs. The subprograms developed here include only the longitud- 
inal motion; however, any number of other subprograms may be added with- 
out revision. The program language is defined in Table VI. All input data 
must be in the foot-pound-radian system unless otherwise noted. Standard 
programing proceducess for FORTRAN 63 are used, [l] The subprograms and 
a sample print out are shown in Figure 35. 

Several of the subprograms require the use of another subprogram to 
compute the roots of the generated polynomials. These subroutines are of 
a standard type but the usage will vary depending on the installation. The 
subroutine used in these programs to compute the roots is RTPLSUB which 
is a library routine available at the Computer Facility of the U. S. Naval 
Postgraduate School. However, any subroutine designed to compute -roots 
of polynomials could be used by rewritting the call statements. 

7.1.1 Main Program - AEROFUN 

The main program may be of any desired format so as to utilize the 
necessary subprograms and print the results. The main program used here 
is dimensioned to take up to 50 sets of d&ta for different flight conditions. 
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The variable NN is used to indicate the number of sets of data for analysis. 

7.1.2 Subprogram COMPTDER 

This subprogram computes those aerodynamic derivatives which must 
normally be approximated from the wind tunnel data. The equations utilized 
for this program are listed in Table VII below. [5] 

TABLE VII 

EQUATIONS UTILIZED IN COMPUTER PROGRAM COMPTDER 



be 

b< 

Cl-qw 

Cy*<j W = — C LC j w |x cul — 



3/ 2 . 

£7 COS A C 



(.sas + e- 385Cht/c) 

A t + 2 (x at “ , X C g)C Lat C L<?t = 2. 5 w ~c~Vt C 



Lot t 



-c 



/R + 2cosA , 


/R 2 TflW 8 A 


1 /K + 2 cos A I" 1 


8/1? + 


2* 


\ A? + 6 cos A 1 



L<* 

= J? t 

^**1* = c. 



Se 1 


Ue\ 


Sw 


a 



[/o] 



C*, 6 =-.o* 



Oy«<j Cyv\<j^ "P 



c. 



Cu s .= 



r 



- it^ T c ^ t 

„ 5. 



S w 






YD 



Se 



Al_ Jh.71 T c 
S w c /*' 



c«u- lr(fl + BC L * t ) 
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INPUTS. , C^, S W' x ac* x cg* 

b€ 

^ £>< • 

NOTE: If Cl^ is based on the wing reference area, assign S w and S t values 
of one. 

be 



b«- 



may either be an input or it may be computed in the program. If 



it is to be an input, assign the value to DEDAI. If it is to be corner 
puted, assign zero to DEDAI. 

OUTPUT: C Lq , C LSe , C mq , C m ^ e , C m * e , C m<< * . 

7.1.3 Subprogram ENGDIMEN 

This subprogram computes the English dimensional derivatives from 
the NACA non dimensional derivatives. The relationships used in this sub- 
program are given in Table II. 

INPUTS: U c , p , W, S w , B, c, fx , z T . l t , C L , C^, C Lu , C L(J , C L4e , 

C D. c d«> C DU' c Dq- c Ds«- C m*' C m* ' C mq' C ms.' C T' 6 o- 
OUTPUTS: X u , X w , X q , X^ , Xp, Z u , Z w , Z q , Z^ , Zip, M u , M w , M^, M q , 

M^,, M t , g cos0 o , g sin 0 O . 



7.1.4 Subprogram STAB 

This subroutine computes the coefficients of the characteristic equa- 
tion, the roots of the equation and the characteristics of the short period 
and long period motions. The equations used in this program are given in 
7.2. 

INPUTS: U Q , 0 Q , X^, X^,, Z u , Z w , Z q , M u , M w , M^, M q . 
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OUTPUTS: (as defined in Table VI) U, V, CONV, WNLP, WNSP, ZLP*, 



ZSP, TLP, THLP„ THSP, TTLP„ TTSP. 

NOTES: Requires a subprogram to evaluate the roots of the generated poly- 
nominal. The coefficients of polynominal are stored in the dimen- 
o. sioned variable C. The order of the polynomial is assigned to the 

variable N. The outputs of the called subroutine (the roots of the 
polynominal) are assigned to U and V and are therefore preserved 
for subsequent calculations. 

7.1.5 Subprogram FORCEFUN 

This subprogram computes and prints the coefficients and the roots 
of the polynominal numerators of the transfer functions c</Tj, u/T^/OC/T, 
u/T. The equations utilized in the program are given in Z. 2. 

INPUTS: U Qf 0Q ff X u , X w , X^ 1 , Z Uf Z w< Zq, Zyy , M u , M w , Mq, 

My, . 

OUTPUTS: None. 

NOTES: Requires a subprogram to evaluate the roots of the generated poly- 
nomial. The operation is the same as for STAB except the output 
is printed internally. 

7. 2 Derivation of Airframe Transfer Functions 
The airframe transfer functions are derived from the equations of 
motion written in Laplace Matrix form given previously as Eqs 2-6. The 
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TABLE VI 





DEFINITION OP SYMBOLS 


Program 

Mnenpnie 


Parameter 


AR 

i 


FR Aspect Ratio 


B 


B (Ijy, SlUg-ft 2 ) 


C, CH 


c chord 


CD 


C D 


CDA 


C D* 


CDDS 


°D*> 


CDQ 


S 


CDU 


C Du 


CL 


C L 


CLA 


C L* 


CLAT 


C L c< t 


CLDE 


C LS 0 


CijQ 


C Lq 


CLU 


c Lu 


CMA 




CMAD 




CMDD 


3° 

Ok, 

CD 


CMDE 


^m5e 


CMQ 


c mq 


CT 


c T 


COOT 


Convergence indicator for roots 
output of STAB only 


D 


Zm Perpendicular distance between 
the* eg and thrust axis. Pos- 
itive for thrust axis above eg 


' 


7F 



TABLE VI(COHT) 



DEDA,DEDAI 


2)g- 


ET 


effectiveness of horizontal tail 


GC THETA 


g cos© 


GS THETA 


g sin© 


HT 


h^. height of root chord position 
of horizontal tail above wing 
root chord. 


m 


Number of sets of data - Main 
Program only 


R 


p (slugs/ft^) 


ST 


s t 


sw,s 


S v; or reference area 


TATJ 


yf- 

1 DSc 


THETAD 


@o ( d °s) 


TLP 


Period of Phugoid Motion 


TSP 


Period of Short Period Motion 


TELP 


Phugoid time to •§■ 


THSP 


Short Period tim to -| 


TTLP 


Phugoid time to .1 


TTSP 


Short Period time to *1 


TR 


taper ratio 


u,uo 


U 0 reference or steady state 
airspeed 


U 


Out out of real parts of roots 
STAB only 


V 


Output of imaginary parts of 
roots - STAB only 


w 


weight 


WNLP 


Wn phugoid 


WNSP 


w n short period 



% 






/ 

TABLE VI (CONT) 



XETA 






XLAMBDA 


j\. sweep (deg) 


XL' 


It 


distance from eg to aero- 
dynamic center of horizontal 
tail 


XMU 


»u 




XMETA 




-\ 


XMQ 


«q 




XMT 






XMW 


% 




XMWD 


% 




XU 


Xu 




XQ 






XT 


x T 




XIV 


*w 




ZE 


angle between thrust line and 
body axis, positive if upward 
(deg) 


ZET 






Z Q, 


Z q 




ZT 


Zy 




ZU 






zv; 


^7 




ZLP 


f of phugoid 


ZSP 


f of 


short period 


A 2c B 

/ 


COMPTDER only - values from 
Perkins & Hage, Pig 10-10* 


XAC 


x a c 


position of aerodynamic center 
of wing in $ MAC 


XCG 


x cg 


position of center of gravity 
in % MAC 
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characteristic equation is the expansion of the main determinant as shown 



below. 



S -X, 



CH = 



a 



- U 0 Xv 



S -Z v 



u« 



g cos 0c 
+ l\ S + 



Z ? , i\ c > s 51 N Q c 



Uo 






S‘ 



M,S 



C H = S’ + A S 3 + 8 S ! + CS + D 



A = -X u -Z w -M, -M*(Z, + U) 

B = X„Z W + M,( Xu+ Z w ) - X« Z„ - Mw( Z, + Uj + Z, + U.) 

+ g sin 0 O 

C = M<j(“X u Z w + X w zj + M u [g co s 0 O _ X w ( U 0 + Z^)] 

+ M^jZ^g cos© 0 - X^g sin© 0 ) + fi v [X u (Z^ +- Uj + g sin 0 o ] 

D = g cos0 o (M w Z a - Z^Mu) -+- g sin© 0 (M u X w — Xal^w] 

The transfer function relating airspeed to a thrust input perturbatidm 
is thens 

0 

_u_ Z 

T C H 



-u 0 x v 



g cosGo 



S - Z w 

-u 0 (m w - SMj 



U 



2 - + 1 



S + 



^ S I N 0 Q 
Uo 



s - n^s 
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LL 

T 



X T (S 3 + A ut s 2 + Bu. t $ + Cu T ) 

C H 

Aut = — ~ Z w — 4* U 0 ) 

But ~ ~ McjZ* — M w ( ^9 + Uo) + si n0 o 

+ -^-[x w (z, + u 4 - 3 c ° s e 0 

m r 

Cut “ M w J5Jn0 o + *^'1 Z w ^ COS 0 O — X w ^ SIN 0 O 



The other transfer functions follow in a like manner. 



oC _ 
T ~ 



CO 

1 

X 

£ 


X r 


^C0S© o 


Zu 

Uo 


0 


_ (uT +1 ) s + 


q SIN Qo 
Uo 


Mu 


n T 


s 2 - 


n,s 



CH 



cC A t S **" BtS + Cj 

~~ Uo(CH| 



= k(s' + C,S + c 2 



K = 



Uo 



. C, = 



By 

At 



C* = 
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>1 O 



a t = x T z u + m t (z^ + U 0 ) 
b t «[riu(z,+ a) - z*mJx t 

-Mt[x ul (Z^ + u 0 ) + 9 -S i n 6q 
C-r = g SIN 0o(Xu M t - fiuX T ) - Z a M T g cos 0 O 





0 


-U.Xv- 


<g COS 0 O 




_Z*L 

U, 


s-z„ 


[Z? +1 ] s+ gsiwe* 

l Uo 1 Uo 


a _ 




-U(Mw + S M*) 


S 2 4- Mq S 






CH 



a 

V 



Zy, X w ( S 2 + A u ^ S 4- But;) 
C H 



Au'r, 






Zq 4- U c 



g cos 9 0 
Z^X W 



My, + M*Z^ 



+ Me 



B 



«n 



^ COS 0c 

~z 7 >c 7 



n n z w - rLz„) + -=X~ y s i n e c 
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QC 
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u 0 
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(j cos 0 O 
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u e 



CH 



Uo 

s 2 -r^s 



Zyjs 3 + A^S a + B^S + 
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Hi 






Z?+ U c 
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7. 3 Evaluation of the Closed Loop Transfer Function. 

The closed loop transfer function was derived as Eg. (2-8) 

T c T f T u 

u = °c T c Tf T 

u c 1 + Dtt lLlLl_ u_ . 

u u e T c Tf T If; Tf T 

If all delays are assumed to negligible, this reduces to 

JEc .Si 

u = a Cj T 

u c l + h _ Is fS 

u 6 T oce T 

Substituting the transfer functions defined in Table IV with 

u _ UN _ *N 

T CH T CH 

gives 

K U (S + UN 

u _ S CH 

u c 1 ."^(s + kn) UN ia 

1 S CH Kx CH 



therefore 

u_ = K U (S + k u )(UN) 

u G S(CH) + K u (S + k u )(UN) - K* (<x:N) 

The denominator of this expression is a 5th order polynomial shown below. 

+ C^S^ + C 2 S 3 + CgS^ + C^S + C 5 where, using 

the same notation as 7 . 2 



C 1 = A + K'yX T 

C 2 = B + K u X,A uT + K u k u X T - K* A*j 

Cg = C + + K u k u X-pA. u >j' - Koc Aoc T 
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C 4 - D + K u X T C uT + K u k u Xi«B u <j- - K otA ^ jC ^ y 
C 5 = K u k i^T C uT 

Then 

u _ XxK u (S + k u )(s3 + A u tS 2 + jl u xS + C u x) 

^ + Cj s4 + C 2 s3~+ C 3 s2 + C 4 S + C 5 

From this expression the root locus was obtained by digital computer means. 
The integral gain produces one variable zero. The other zeroes of the func- 
tion are invariant with changes of the gains. The root locus plots obtained 
are shown in Figures 33 and 34. The most desirables gains can now be 
chosen from these plots using standard pole-zero methods. [2] However, 
this type of an analysis would have to be made for each loop involved to .. 
arrive at the best combination of gains for both the airspeed and the angle 
of attack response to grists. 



7 . 4 Analogue Simulation 

Substituting the aerodynamic values into Eqs 2-4 gives 

u = -. 060U <*- 3. 32 - 31.850 +.00145T - 1.647? 

* = -. 00113u - .4265« + 0 - . 019350 - .08157? 

0 = 1. 85x10“ 4 u - 1.14oC _ .0415* - .3380 - 4. 55 x10"6t 
- 2.2577 

These equations were then scaled using the relations defined in Table 
VIII. Real time was used. Using this scaling^ the voltage form of the equa- 
tions becomes 

u = -. 101u - . 0116 ^ - . 1110 + . 292T - . 00677 
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°C = 544u - •4265 oC + 0 - . 019359 - .08157? 

f = . 0311u - „ 398 - .0145^- . 11810 - .7857? 



09 IT 



The power compensator equation was scaled as (Eq 2-5 without 
delays) 



The delays were simulated by the general relationship shown below 
were 0 O is the delayed signal and 0^ is the input signal. 



The values of 1/T can be obtained by either potentiomenter settings or by 
adjusting the input gains. 

The schematic diagram for the wiring circuits is given in Figures 7 
and 8. The number on each input is the gain of the, input where the gain is 
as defined below: 

For a summing amplifier For an integrator 




118.9 Ue 




1 + TS 



1 



then 0 O = 1/T0i - 1/T0 O 



or 




which gives the simulation circuit 





C 
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Gain 



Rf 

Ri 



Gain 



RiC 



The potentiometer settings are given in Table IX. In order to keep 
the wiffirigg diagram as simple as possible, feedback circuits are indicated 
by the use of numbered potentiometers or numbered units shown dashed. 

The squares indicate recorder channels. SW 1 and SW 2 are internal double 
throw-double pole switches. NC refers to the "normally closed" position 
of SW 2. Step gusts were simulated by means of initial conditions on the 
airspeed and angle of attack integrators. Sinusoidal gust inputs were 
obtained by means of a frequency oscillator. 
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TABLE VIII 
SCALED VARIABLES 


. (9# 

Jjjr 


Parameter 


Maximum Value 
(per 100 volts) 


Voltage 

Equation 


Scale 


u 


50 ft/sec 2 


d = 2 u 


.5 ft/sec/volt. 


U s u c 


81]. .3 ft/sec 
56 kts 


u = l.l89u 
u = 2u 


•81j.3 ft/sec/volt 
.5 kt s/volt 




.1714-5 rad 
10 deg 


* = 573 *■ 

c*. = 10 <*• 


.00171l 5 rad/volt 
.1 deg/volt 


»• 

© 


.5 rad/sec 2 


¥ = 200 e 


.005 rad/sec 2 /volt 


e,e,s-,\ 


.1714-5 rad 
10 deg 


1= 573 e 
© - 10 © 


.00171 {.5 rad/ volt 
.1 deg/volt 


T 


1000 lbs 


7 - T/100 


100 lb/volt 
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TABLE IX 

POTENTIOMETER SETTINGS 



5t # 


Setting 


Input 

Parameter 


Gain 


1 


*292 


T 


1 


2 


.006 


VI 


1 


3 


.101 


u 


1 


k 


.0116 


cC 


• -1 


5 


.111 


0 


1 


6 


.101/. 0116 


Ug/oCg (Gusts) 


1/1 


8 


.0815 


'a > 


1 


9 


.272 


u 


2 


10 


.2133 


0 c 


2 


12 


.019k 


e 


1 


13 


.272/2133 


uV«<-„ (Gusts) 

f=> s 


2/2 


it 


.091 


T 


1 


15 


.393 


v\ 


2 


16 


.031 


u 


1 


17 


.199 




2 


18 


.Olij.5 




1 


21 


.118 


• 

0 


1 


22 


.031/0 


u -/<£'„ (Gusts) 

^ O 


1 


23 


Various 


u c 


1 


2t 


0/.199 


Ug/cC g (Gusts) 


2 


27 


Eu 


u e 


10 


26 


K 




1 


29 


KA 


f u e 


1 
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TABLE IX ( C Oi'IT ) 



31&32 

V 


.5 


T c (Throttle Servo) 


10 


3 3&3li 


.hi5 


Tf (Engine Delay) 


2 


35 


Various 


X 


1 


36 


Various 


Ac 


1 


38 


.5 


X delay 


1 
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